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00 '. Abstract 

O 

^D I We estimate the single-spin asymmetries (SSA) which provide the access to transver- 

^^ ' sity as well as to Boer-Mulders and Sivers PDFs via investigation of the single-polarized 

^ . Drell-Yan (DY) processes with pp, pD and DD collisions available to RHIC, NICA, COM- 

PASS, and J-PARC. The feasibility of these SSA is studied with the new generator of 



:z; 



^ 



s 



polarized DY events. The performed estimations demonstrate that there exist the such 
kinematical regions where SSA are presumably measurable. The most useful for PDFs 
extraction are the limiting kinematical ranges, where one can neglect the sea PDFs contri- 
butions which occur at large values of Bjorken x. It is of interest that on the contrary to 



Q^' the Sivers PDF, the transversity PDF is presumably accessible only in the especial kine- 



cih: 



matical region. On the contrary to the option with the symmetric collider mode (RHIC, 
5^1 NICA), this is of importance for the COMPASS experiment and the future J-PARC fa- 

r^,' cility where the fixed target mode is available. 



^ . PACS: 13.65.Ni, 13.60.Hb, 13.88.-fe 

O 
00 

^ 1 Introduction 

(N : 

C^ ! In this paper we focus on DY processes in collisions of transversely polarized protons and 

00 ! deutrons, providing us an access to the very important and still poorly known sea and valence 

O ■ transversity, Boer-Mulders and Sivers PDFs in proton. At present the such processes are 

K*" ■ available to RHIC [T]. Their studies are planned at J-PARC [2] facility and, in principle, are 

S^ - possiblqflat COMPASS, where unpolarized proton beam and both polarized proton and deutron 

targets are available. Besides, at present JINR starts the new NICA/MPD project based on 

the development of the existing Nuclotron accelerator for the new facility creation: the heavy 

and light nucleus collider NICA |4], [5j. In particular, the possibility is now considered to study 

the collisions of the polarized proton and deutron beams at the second interaction point (IP) 

at NICA. It is of importance, that all these experiments do not duplicate each other, but are 

complementary, providing us by the the information on the different PDFs measured in the 

differenlo kinematical regions. 



^e-mail address: sisakian@jinr.ru 

^e-mail address: shev@niail.cern.ch 

•^e-mail address: nagajcev@sunse.jinr.ru 

■^e-mail address: ivon@jinr.ru 

^At present the Drell-Yan program at COMPASS focuses on the pion-proton (deutron) collisions (see, for 
example, [3]). However, the possibility to study DY processes with pp and pd collisions is now also under 
discussion. 

^ For pp collisions the center of mass energy is 200 GeV for RHIC. For COMPASS it can be varied from 20 
to 27 GeV (upper bound corresponds to 400 GeV primary proton SPS beam). For J-PARC y/s is planned to be 
about 8 GeV at first stage and about 10 GeV for the second stage. For NICA it is planned to be about 20-26 
GeV. 



The leading twist kT integrated transversity PDF AtQ = hig, as well as the leading twist 
unpolarized q = fiq and longitudinally polarized (helicity) Ag = g^ PDFs, is of the crucial 
importance for understanding of the nucleon spin structure (see Ref. [6] for the comprehensive 
review). At the same time, nowadays the study of quark transverse momentum kx depen- 
dent PDFs is also among the special issues in hadron physics. Of particular interest, are 
two leading-twist T-odd kx dependent PDFs: Sivers function /^/(x, /c|.) and Boer- Mulders 
function h^(x,k^). While Sivers function represents the unpolarized parton distribution in 
a transversely polarized hadron, the Boer-Mulders function denotes the parton transversity 
distribution in the unpolarized hadron. 

At present the Boer- Mulders PDF is still not measured, while the Sivers [3, E] and transver- 
sity ^ PDFs were (preliminarily, with rather big uncertainties) extracted from the SIDIS data 
collected by HERMES ^lOj and COMPASS ^Aj collaborations. At the same time the analysis 
of SIDIS data suffers from the poor knowledge of the fragmentation functions, and especially 
it concerns the Collins fragmentation function which is necessary for the transversity extrac- 
tion [H]. In this respect the Drell-Yan processes possess the essential advantage since they 
are free of any fragmentation functions. Besides, DY measurements should accompany SIDIS 
measurements to check the important QCD prediction [12] (see also ^iBj and references therein) 
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for T-odd PDFs j^t ^^"^ ^i"- ■'■^ i^ relevant to notice in this connection that while the Sivers 
PDF was already extracted from SIDIS (with rather bad precision but at least the sign of 
ffV is seen) there still was not the respective analysis for h\- . Recently in the 

•' ^-' SIDIS ' i- J I SIDIS 

paper [14j the estimations were presented of the possibility to extract the azimuthal asymmetry 
(cos 20) (giving the access to h\^ from the combined analysis of all existing and planned SIDIS 
measurements. Thus, if this SIDIS program would be realized, then DY measurements can 
allow to check sign change effect for both Sivers and Boer-Mulders PDFs. 

It is well known that the double transversely polarized DY process H{Hl -^ l^l^X allows 
to directly extract the transversity distributions (see Ref. [6] for review). However, at the same 
number of collected DY events, the double spin asymmetries suffer from the much more large 
statistical errors (product of the beam and target /beam polarizations in the denominator) than 
the single spin asymmetries. This is especially important for the double-polarized DY processes 
with pp (as well as with pD and DD) collisions, where because of the small value of transversity 
PDF for the sea antiquark in proton (neutron), the double spin asymmetry is estimated to be 
a few percents maximum [TB]. Thus, in the case of pp, pD and DD collisions we focus on 
here, the double polarized DY processes seem to be not too useful and we need an alternative 
possibility allowing to extract the transversity PDF from the combined analysis of unpolarized 

H^H2 -^ l^-l- X, (2) 

and single-polarized 

H^hI -^ l+l- X (3) 

DY processes. Besides, namely unpolarized and single-polarized DY processes give us also an 
access to Boers-Mulders and Sivers PDFs, which are very intriguing and interesting objects 
in themselves. On the other hand, in the processes ([2D and (^ the access to PDFs we are 
interesting in is rather difficult since they enter the respective cross-sections [TB| in the complex 



convolution with each other, so that at first sight it is impossible to avoid some models on 
the kr dependence of PDFs. To solve this problem the qt weighting approach [TTl [HI [191 
[20] was recently applied in Ref. [7] to Sivers effect in the single-polarized DY processes ([3j), 
and in Refs. fHX [22] with respect to transversity and Boer-Mulders PDF extraction from 
both unpolarized and single-polarized DY processes ([2]) and ([3]). In two last cited papers 
we considered the DY processes with antiproton-proton and pion-proton collisions. At the 
same time the DY processes with the proton (deutron)-proton(deutron) collisions are also very 
important since they provide the access to sea PDFs. Within this paper we will estimate both 
types of single-spin asymmetries (SSA), which give us respectively access to Sivers PDF |7l [8] 
and to transversity and Boer- Mulders PDFs [211 ES] . At first sight it seems that DY processes 
with proton(deutron)-proton(deutron) collisions are strongly suppressed because there is no 
valence antiquark in the initial state there. However, we will see that there exist the kinematical 
regions where both SSA take quite considerable values. 

2 Transversity and T-odd PDFs via Drell-Yan processes 
with pp collisions 

The procedure proposed in Refs. [211 122] allows us to extract from the processes ([2]) and ([3]) 
the transversity hi and the first moment 

hi^'\x) ^ J d'kT (^) hi^{x„ k|) (4) 

of Boer-Mulders hi PDF directly, without any model assumptions about /cr-dependence of 
hi{x^k^). Applied to unpolarized DY process ([2]) with pp collisions this general procedure 



giver' 



pp^m-x E,e2[/i,(xi)/i,(x2) + (g^g)] ' 



(5) 



where k is the coefficient at cos 20 dependent part of the properly qt weighted ratio of unpo- 
larized cross-sections: 

^- T^v^w ' ^6) 

^= -|-(7(l + cos2^) + fccos20sin2^). (7) 

At the same time, in the case of single-polarized DY process ([3|), operating just as in Ref. 
one gets 

1 _ 1 Eqel[htg^\xp)hiq{Xpi) + {q^ g)] 
2 Egel[fiq{xp)fiq{Xpi) + {q-^q)] 



^Eq. ([5|) is obtained within the quark parton model. It is of importance that the large values of coefficient 
at cos2(/) in the ratio of DY cross-sections can not be explained by the leading and next-to-leading order 
perturbative QCD corrections as well as by the high twists effects (see [TB] and references therein). 



where the single spin asymmetry (SSA) Ah is defined as 

, _ Idnd<PsJd\TihT\/M,)smi^ + <Ps,)[daiS2T)-da{-S2T)] , . 

' J dndcPsJ d''ciT[dcr{S2T) + dai-S2T)] ' ^^ 

All the notations used are the same as in Ref. [21] (see Ref. [6] for details on kinematics in the 
CoUins-Soper frame we deal with). 

sin{<f)—<ps) tif 

Notice that SSA Ah is analogous to asymmetry Aijrp ^ (weighted with sin(0 — 0s) and 

the same weight qt/Mn) applied in Ref. [7] with respect to the Sivers effect investigation in 
the single-polarized DY processes. For DY process pp^ — ^ l^l~X we study here the expressions 

sin((A— (Ti<;)-rr— 

for A^^ "''^ look as (see Eqs. (14), (15) in Ref. ) 



sin(0-<^,)^ ^ / dndct>s, I d^TJM/M,) sin(0 - C^s,)[dcT{S2T) ' dajS^T)] . . 

"""^ U dnd^sj d^qT[da{S2T) + da{-S2T)] ' ^ ^ 



,--('^-</'s)iS^ _ ^ E,e^[/iT^'^'(v)/ig(a:p) + (g ^ q) ] 
Eqel[fiq{Xp^)fig{xp) + {q-^ q)] 



A^T ^ = 2 ^1 ''':ti / '::'"" , — r^, (111 
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where 

ft^"\x) ^ J d'^kr (^^ J f^r'i^, k^) (12) 

is the first moment of the Sivers function fij!{x, k^)- Notice that factor 1/2 in denominator of 
Eq. (TTUj) (see also Eq. (7) in Ref. [23] ) was introduced in Ref. [7] (where the Sivers effect was 
studied in both SIDIS and DY processes) for consistence with the respective semi-inclusive SSA 
studied by the HERMES - see ^J and references therein. Since within this paper we also will 
study SSA given by Eqs. (fTOj) and (ITTl) . for comparison purposes it is convenient to introducqfl, 
by analogy, SSA 

^^ Eqel[fiq{xp)fiq{Xpr) + {q^q)] 



Let us now consider SSA given by Eqs. (ITTll and (1131) . On the contrary to valence PDFs, 
the sea PDFs dominate at small x and rapidly die out when x increases. That is why in the 
case of pp^ collisions we deal with the regions are of importance where the Bjorken x for sea 
PDFs take the small values, while, by virtue of the relation 

XpXpi; = Q'^/s, (14) 

the valence PDFs occur at large x values. Indeed, in such the regions we can neglect the 
contributions to SSA containing sea PDFs at large x (later we will see that this is really good 
approximation) and, thereby, to essentially cancel the number of extra unknown PDFs entering 
the asymmetries. 



^ Certainly, from the practical point of view, the such rescaling is not especially useful: when the asymmetry 
is multiplied by a number the error is multiplied by the same number too. However, it is convenient to consider 
both studied SSA given in the same scale in order to estimate and compare their feasibility at the same statistics 
of DY events. 



Thus, let us consider two limiting cases Xp ^ Xpr and Xp <^ Xpt . In the first case 
neglectingjthe terms containing the sea PDFs at large Xp, one arrives at the simplified equations 
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Xp^PX^^ 



4/lT_ {Xp^)flu{Xp) + fiT {Xp^)fld{Xp) 
4/ln(v)/l«(^p) + hd{Xp^) fid{Xp) 



(16) 



and 
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sm{<p+<t>s)^ 
UT 
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Xp>Xp| 



4/iiJ ' {xp)hiu{Xpi) + h^} ixp)hld{xp^) 
4/i«(a;p)/i„(xpT) + fid{xp)fid{xp'\) 



(17) 



Then, taking into account the quark charges and u quark dominance at large x, Eqs. ( TT6l) and 
(TT7|) are essentially given by 
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and 
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hiu {xp)hlu{Xp^) 

'^p»V ~ fluixp)fluixp^) 

Analogously, in the second limiting case 

Xunpol "^ XpQi, 

one gets 
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^, 2 ^fiT {xp^)fiu{xp) + fiT {xpr)fid{xp 
^P«V 4/i„(xpT)/i«(xp) + fldixp^)fldixp) 



and 
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4/iiJ n^p)^i^t(a^pT) + /^li ixp)hld{xp^) 
4^fiu{xp)fiuixp^) + /id(a;p)/id(xpT) 



with d quark contributions, while 
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(19) 
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(24) 



^Notice that all over the paper we neglect in our calculations the contributions of strange PDFs which 
produce really tiny corrections. 



neglecting d quark contribution. 

Later (see Section 3) we will see that even the double approximations given by Eqs. (TT9!) and 
(!24|) are in a very good agreement with Eqs. (TTTl) and (TT3|) . It is of importance since allows us 
to cancel out the extra unknown variables entering the equations for measured asymmetries. In 
particular, this give us the interesting possibility to extract the ratios hiu/h^!^ and hiu/hi^ 
directly, without application of the fitting procedure with a set of assumptions on extra unknown 
variables. Indeed, let us return to the unpolarized DY process with pp collisions, Eqs. ©-([Tj). 
In the limiting cases Xi ^ X2 and Xi <^ X2 Eq. ([5]) is reduced to the equations 
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and 
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^l»^'2 fluiXi)fiu{x2) 

^,o hi^^\xi)hii^\x2) 

^l«^2 flu{Xi)fiu{x2) 



(25) 



(26) 
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respectively. Then, having in our disposal the quantities k{xi,X2) and Ajjrp 

Xi,Xpt = X2) measured in unpolarized and single-polarized DY processes, and, combining Eqs. 

( fT9|) . ( l24l) with Eqs. ( !25l) . ( l26l) . we can obtain the ratios hiu/hi^ and hiu/hiu from the 
equations 



^sin(ih+<l>s)4f- ^sm(6+(hs)Tr-- 



^UT 



hti^\xi) k -i»-2' ht!:^\xi) k 

3 Estimations on SSA in pp collisions 



xi<SiX2 



(27) 



sin{(t)—<f>g) ?J 

Let us first estimate SSA A^r^ ^ 



given by Eq. ([n]) . Notice that for the RHIC 



kinematics SSA Jt™jC' ' weighted only with the projecting factor sin(0 — (\)s) was already 
in detail studied in Ref. [23]. This asymmetry within the Gaussian model applied in [23] 

sin(0— </is)-j^IL 

is just proportional to the qr weighted SSA Ayj, ^ which we consider in this paper: 

^sin{</,-</,s) _ ^gy ^^ A^!f~'*"'^^ {a%l^ss - 0-81 " (1 ± 10%)), SO that we need not to repeat here 
the calculations for the RHIC kinematical domain. Instead, we present the estimations for the 
NICA kinematics, where two colliding 10 ~ ViGeV proton beams assumed to be available [1], 
[5]. We perform the calculations for Q^ values below and above J/\E' threshold Q^ = 9.5 GeV^. 
For the estimations we use three different fits for the Sivers function: fits I and II from Ref. 
[7] and also the latest fit from Ref. [5], which we denote as fit III. These fits are given by the 
following parametrizations 

FitL xf^^^^"" = -xf^^^^'' = OAx{l-xf, (28) 

Fit II: xf^r^'^"" = -xftr^^^" = 0.1x°-3(l - x)\ (29) 

Fit III: xft^^^"" = -x/ir^'^'' = (0.17...0.18)x°-''^(l - xf. (30) 

For the first moments (TT2!) of the sea Sivers PDFs entering Eq. (TTTl) we use the model (with 
the positive sign) proposed in Ref. ^23j (see Eqs. (10) and (11) in Ref. [25]): 

ft^'^\x) h^{x) + h,{x) 



ft^'^\x) fluix) + fUx) 



(31) 



For the unpolarized PDFs entering Eq. (fTTj) we use GRV94 [24j paranietrization. The results 
of estimations for the different Q^ values are presented in Fig. [H 

Looking at Fig. [T]one can see that the asymmetry takes the largest values near zero value 
of Xp — XpT and when this difference becomes positive. 

Notice that besides of parametrizations I, II and III on the Sivers PDF there exist also 
the parametrizations from Refs. ^25j, [26j. Our estimations^ with these parametrizations 
show that the values of A^^^ are quite similar to the respective values obtained with the 

parametrizations I, II, III in the region Xp > Xpt and even higher in the region Xp < Xp\. 

Thus, the prediction on A^^ with parametrizations I, II, III could be considered as even 

underestimated predictions: if Jt^j; ' predicted by parametrizations I, II, III will be seen 
within the errors, then the larger SSA predicted by the fits from Refs. |25j and [2S] will be 
definitely measurable too. 

Our calculations performed for the COMPASS and J-PARC kinematics (in comparison with 

sin((i!)— (/)g) ■ ,^ 

RHIC rather close to the NICA one) produce for SSA Auj. ^ the plots very similar to the 

sin((/i— 0s)-j^^ 

respective plots for the NICA kinematical range. The such behavior of A^jj. ^ encourage 

one that it can be measured not only in the collider mode (RHIC, NICA) but also in the fixed 
target experiments (COMPASS, J-PARC). The point is that, as a rule, in the experiments with 
a fixed target the acceptance of the detector allows us to register mainly the events with the 
positive Xp: 

Xp ^ Xij(,am ^target ~ U, \'^^) 

where Xbeam and Xtarget are the Bjorken x values for the quarks inside of the beam and target 
protons, respectively. Thus, the option Xp > Xpt can be realized by both COMPASS and J- 
PARC facilities, where unpolarized proton beam and polarized proton target can be available. 

sin(0— (/15) ■ ,^ 

On the hand, the region Xp < Xpt, where SSA A^rp ^ is also quite considerable, definitely 

can be reached by the collider experiments RHIC and NICA, and, presumably, by J-PARC, 
where the option with the polarized beam is also planned (see Ref. 0). 

Let us now estimate SSA A^j, ^ given by Eq. (IT3l) . Since neither the Boer- Mulders 

function nor its first moment are still not measured, we will use in our calculation the Boer's 
model (Eq. (50) in Ref. [H]) which produces the good fit for the NAlO [2B] and E615 [2S] data 
on (anomalously large) cos (20) dependence of DY cross-sections. This model gives for the first 
moment (jl]) entering Eq. flT3l) the value hiq (x) ~ 0.163/i(x). We also apply the following 
assumption for the first moment of the sea Boer- Mulders PDF 

ht^\x) UxY ^''^ 

Notice that this assumptions is similar to the assumption fl3T|) . however, namely Eq. fl33|) is 
consistent with the Boer's model which we use here. 

Recently, for the fist time, the transversity PDF was extracted [H] from the combined data of 
HERMES pro], COMPASS [30] and BELLE JET] collaborations. However, because of the rather 



^"^Dealing with the paranietrization from Ref. [25j one should remember that the notations for Sivers PDF 
^^ Iq/H^ and /j^ differ by the sign and extra multipher: f^{x, k|,) = — (A//2|kT|)A^/q /^t {x, k|,) (see "Trento 
conventions" in Ref. |i27j ) . 



poor quality of data the error band surrounding the fit on hi is very large (see Fig. [2]), and, 
besides, the authors of Ref. [9] were compelled to apply the large number of approximations. 
In particular, the approximation of zero sea transversity PDF was applied. However, as it was 
stressed before, in the case of proton-proton collisions namely the sea PDFs play the crucial 
role. That is why here we will apply two versions of evolution model for transversity instead of 
the fit from Ref. [9]. First is the model where the Soffer inequality is saturated [15j: 



1 



2 



1 



hi,ix,Qi) = - g(x,Q^) + Ag(x,Q^) , hig{x,Qf,) = - g(x, Q^) + Ag(x, Q^) (34) 



2 



at low initial scale (Qo = 0.23GeV^), and then hiq, hiq are evolved with DGLAP. Certainly, 
this model rather gives the upper bound (maximal value) on SSA. In the second version of 
evolution model (see [6l |32] and references therein) the valence and sea transversity PDFs are 
assumed to be equal to helicity PDF Ag at the same initial scale (model scale) and then hiq and 
hiq are again evolved with DGLAP to the required Q^ values. Certainly, this model is much 
more realistic one because at the model initial scale a lot of models predict [B] that hi = Aq. 
It is of importance that the curve corresponding to this version of evolution model lies just 
inside the error band for the fit of Ref. [9] - see Fig. [2l Thus, this version of evolution model 
is consistent with the analysis of Ref. [9] . 

We present here the estimations of SSA A^J' ^ below and above J/ip resonance for 

the strongly different RHIC (Fig. [3]) and NICA (Fig. H]) kinematical conditions. The respective 
plots for the COMPASS and J-PARC kinematics are again very similar to that for the NICA 
kinematical range. 

Looking at Fig. [3]and|l]one can see that for both RHIC and NICA (as well as for COMPASS 

and J-PARC) kinematics the asymmetry A^j. ^ is negligible at Xp > Xpt and is quite 

considerable at Xp < Xp^. In this second case SSA Ajjrp ^ takes its maximal values (about 

5-10%) when Xp — Xp] takes the large negative values. Thus, one can again see the advantage 
of the symmetric collider mode (RHIC, NICA) where the cases Xp < Xpt and Xp > Xpt do 
not differ especially. On the contrary, for the fixed target mode these cases essentially differ 

1 IT 



sin(<^+0s)i^ 



because of the acceptance restriction fl32l) . Thus, to obtain nonzero A^j, ^ with the 

fixed target, one should either manage to avoid the restriction (132|) dealing with unpolarized 
beam/polarized target option (forward-backward geometry spectrometer), or, in the case of the 
forward-geometry spectrometer (acceptance restriction ( l32l) holds), one should deal with the 
polarized beam/unpolarized target option. Regretfully, both these options are hardly possible 
for the running COMPASS experiment. At the same time, the option with the polarized proton 
beam is now planned at J-PARC facility [2]. 

In conclusion of this section, let us estimate how good are the approximations given by 
Eqs. (IT6l)- (iT9l) in the case Xp 3> XpT, and by Eqs. (l2Tl)-(l24l) in the case Xp <C Xpt. This is 
very important for the analysis since these approximations allow us to cancel out the extra un- 
known variables entering the equations for measured asymmetries. The respective calculations 

sin(4'—<t>s)jf— JTTl __ 

of A^rp '^ are presentecO by the Table [H From this table it is seen that the approxi- 

mations fITB]) and fl^T]) obtained by neglecting the sea PDFs contributions at high x work very 
well (column "B"). At the same time, the agreement of the results obtained with the double 
approximations (fT8|) and (123|) (column "C" ) with the result of column "A" is worse (but still 



^^For better readability of the paper we present here only two tables for NICA kinematics. Our calculations 
show that for all another considered within the paper kinematical conditions the results are absolutely analogous. 

8 



sin{<f>—<ps) M 

Table 1: NICA kinematics. Values of Auj, ^ calculated by using Eq. (Uljl (labeled as 

A) in comparison with two approximations: Eqs. (fT6|) . (12T|) . labeled as B, and Eqs. (TTSl) . (!23ll . 
labeled as C For /u^ the parametrization from Ref. [8j is used. 



Xp — ^T>T yV O Ly Xp — XpT xV 13 L^ 



-0.4000 0.0189 0.0184 0.0277 -0.4000 0.0178 0.0170 0.0277 

-0.5000 0.0131 0.0129 0.0190 -0.5000 0.0132 0.0129 0.0204 

-0.6000 0.0087 0.0086 0.0125 -0.6000 0.0093 0.0093 0.0142 

-0.7000 0.0053 0.0053 0.0076 -0.7000 0.0061 0.0061 0.0091 

-0.8000 0.0028 0.0028 0.0040 -0.8000 0.0033 0.0033 0.0049 



0.4000 


0.0514 


0.0525 


0.0614 


0.4000 


0.0828 


0.0849 


0.0984 


0.5000 


0.0486 


0.0491 


0.0556 


0.5000 


0.0811 


0.0820 


0.0922 


0.6000 


0.0460 


0.0462 


0.0509 


0.6000 


0.0788 


0.0792 


0.0867 


0.7000 


0.0437 


0.0438 


0.0471 


0.7000 


0.0764 


0.0765 


0.0818 


0.8000 


0.0417 


0.0418 


0.0439 


0.8000 


0.0742 


0.0742 


0.0775 



not so bad). This is not surprising since in the applied parametrizations (!28l) -( l30l) for the Sivers 
PDFs [7j, [8] the relation f^f = —fix is used, which is argued within the 1/Nc expansion 
approach [33l [3l] . 

On the other hand, in the case of asymmetry A^j, ^ both kinds of approximations 

([HD, ([IHD and ([22]), ([21]) work very well (see Table [2]). 

Thus, dealing with the kinematical region Xp ^ Xp^ one can safely use the approximation 
( TT6]) and, thereby, get the access to the first moments of the sea Sivers PDFs. Dealing with 
the approximation (ITS]) one should be more careful - if indeed f^ji ~ ~fiT then double 
approximation (I18p is suitable only for some rather rough estimations. At the same time, we 
have no access to the transversity and Boer-Mulders function in the region Xp ^ Xpt since 

sin(0+(/,s)^ 

the asymmetry A^rp '^ is negligible within this kinematical region. Let us recall that for 

the fixed target mode, this kinematical region corresponds to the option with the unpolarized 
beam and polarized target (if, certainly, we have in our disposal only the forward-geometry 
spectrometer - see Eq. ([32]) ). 

Let us now consider another limiting case Xp -C x^t , which in the case of the fixed target 
mode and forward-geometry spectrometer corresponds to the option with the polarized beam 
and unpolarized target. In this case the situation is absolutely opposite. Here the asymmetry 

sm{<l>+(f>s)jf- 

Ajjrp ^ is quite considerable (and presumably is measurable as it was discussed above). 

As it has just shown, in this limit we can safely apply even double approximation ([2^ . This 
give us the interesting possibihty (see Section 2) to extract the ratio hiu/h-i^ directly, without 
application of fitting procedure with a set of assumptions on extra unknown variables. 

4 SSA in pD and DD collisions 

As usual, the inclusion of the deutron beam/target can allow us to find PDFs of u and d quark, 
in separation. 

Applying SUf(2) symmetry to the results of Section 3 one immediately gets the respective 



Table 2: NICA kinematics. Values of Ajjj, ^ calculated by using Eq. ([8]) (labeled as A) 

in comparison with two approximations: (fT7|) . (!22l) . labeled as B, and (TT9|) . (!24l) . labeled as C. 
The evolution model for transversity with /iig(g) = Ag(Ag) at Ql = 0.23GeV^ is used. 



s = AOOGeV^, Q^ =4:GeV^ 
Xp — Xpt A B C 



-0.4000 
-0.5000 
-0.6000 
-0.7000 
-0.8000 



g = 400 Gey ^ Q^ = 15GeV^ 
Xp — Xpt ABC 



-0.0761 
-0.0838 
-0.0894 
-0.0940 
-0.0987 



-0.0800 
-0.0856 
-0.0902 
-0.0943 
-0.0988 



-0.0912 
-0.0948 
-0.0975 
-0.1000 
-0.1029 



-0.4000 
-0.5000 
-0.6000 
-0.7000 
-0.8000 



-0.0783 
-0.0864 
-0.0926 
-0.0980 
-0.1038 



-0.0833 
-0.0887 
-0.0936 
-0.0984 
-0.1039 



-0.0951 
-0.0983 
-0.1012 
-0.1041 
-0.1078 



0.4000 


0.0063 


0.0067 


0.0068 


0.4000 


0.0200 


0.0216 


0.0220 


0.5000 


0.0052 


0.0054 


0.0054 


0.5000 


0.0176 


0.0184 


0.0186 


0.6000 


0.0044 


0.0045 


0.0045 


0.6000 


0.0156 


0.0159 


0.0160 


0.7000 


0.0038 


0.0038 


0.0038 


0.7000 


0.0138 


0.0139 


0.0140 


0.8000 


0.0033 


0.0033 


0.0033 


0.8000 


0.0123 


0.0124 


0.0124 



results on SSA for Drell-Yan processes in pD and DD collisions. For SSA giving an access to 
the Sivers PDF, in the limiting case (fTSjl one gets instead of Eq. (fT6|) the equations 



and 



A 



sin(0-(/>s)^ 



A 



UT 



UT 



(xd > x„u 



Dp^^l+l-X 



^Xp S> X£)t ^ 



pD'i^l+l'X 



while in the limiting case ( 1201) one obtains 



4/i„(XpT) +/ld(XpT) 



(35) 



A 



UT 



,xd >a;BT, 

(i)d. 



DD^^l+l-X 



hu{xD^) + hd{xj:)}] 



(36) 



and 



A 



sm{4,-^s)j^ 



UT 



^UT 



{xd < Xp 



_(!)« 



(l)d. 



Xp -CX^iT, 



Dp^-.i+i-x 4:hu{xpr) + fidixp,) 

sin(</.-0s)^ 



pDT^l+l-X 



A 



UT 



'Mm 



Xd ^ a^DT, 



(37) 



DD'^^l+l-X 



J^T^^'^iXu^) + /i"t^^°'(x^t 



flu{xD^) + fld{xD^) 

At the same time, SSA giving an access to transversity and Boer-Mulders PDFs look as 



A 



sin{<^+0s)^ 
UT 



[Xd > x„u 



M^)i 



AW 



DpT^l+l-X 



[hiu (a^J)) + hd ixD)][4:hlu{xp^) + hid{xp^)] 

[fluixD) + fldixD)][4:fluiXpl) + fldiXp^)] 



(St 



(39) 



A 



UT 



[Xp > XdT, 



-(1)/ 



-(1) 



[4/?-in jxp) + h^d {xp)][hiu{xD^) + hld{xD^)\ 

pD^^l+l-X - [4/i„(Xp) + fld{Xp)][flu{xD^) + fldixDl)] 

10 



(40) 



and 
A 



UT 



{xd > xor 



DDT^l+l'X 

in the limiting case (fT5|) . while 



A, 



UT 



{Xd <^ Xr 



Dp^-*l+l-X 



AW I 



-(1) 



[hiu (xp) + Kd {xD)][hlu{xI:)^) + /j-id(xgT)] 
[fiu{.XD) + /id(a;D)][/i«(a;DT) + /idla;/??)] 



-Ml), 



-Ml) 



[h^^ '{xp) + h^} ixD)][Ahiuixp-i) + hld{xp^)] 
[fiuixD) + /id(a;D)][4/i„(xpT) + fidixp^)] 



(41) 



(42) 



sin(<j>+<f>s) 

Ajjrp 




[Xp 


< 


a^Dt) 


and 










sin((/)+0s) 


9T 


{xd 


< 


Xd^) 



pDl^l+l-X 



[4:htu^\xp) + hfj^\xp)][hiu{xDr) + /^id(a^DT)] 
[4/in(xp) + /id(a;p)][/i«(xDT) + fidixoi)] 



(43) 



-(1), 



;^(i) 



DDt^l+l-X 



[/j,i„ ^(xp) + /ii] {xD)][hiu{xD^) + hld{xD^)] 

IfluixD) + fldixoWluixDr) + fldixD^)] 



(44) 



in the limiting case (120|1 . 

As it was mentioned above (see Section 3) there exist the strong theoretical arguments 
[33] . [3i] based on l/Nc expansion that the sum of the u and d quark Sivers first moments 



f^ and fij! is also very small quantitjo. Besides, the QCD evolution predicts small 
values of the sea transversity distributions even at small x values [32]. Thus, in the case 
of polarized deutron in initial state, almost all respective SSA (see ( 136|) . fl38l) . (HOj) . fHTl) ) 
presumably should be very small quantities (and our calculations confirm it), compatible with 
zero within the errors (certainly, it should be carefully checked by the respective measurements 
at RHIC, NICA, COMPASS and J-PARC). The only SSA which could take considerable value 
are SSA containing the sum hiu^Xor) + hid^Xor) (see Eqs. (H3|) . dBD). The point is that the 
analysis [9] of the COMPASS data [30] obtained on the deutron target produced the possibility 
of nonzero sum hiu + hid- In accordance with this analysis /ii„ and hid are of opposite sign but 
diffeio in their absolute values (see Fig. 7 in Ref. [U]). However, the uncertainties on hiu and 
hid are too large (see the error bands in Fig. 7) to realize is the quantity hiu + hid zero or not. 
Thus, the respective measurements of SSA in DY processes with polarized deutron could shed 
the light on this problem. 

On the contrary to the case of polarized deutron, in the case of polarized proton all SSA 
could take the considerable values. Our calculations show that they are of the same order 
of magnitude as the respective SSA in the case of pp^ collisions - see Figs [5], El Notice that 
for brevity we present the relations between SSA for pp^ and Dp^ collisions only for NICA 
since they are very similar (practically the same) to that for RHIC, COMPASS and J-PARC 



kinematics. In Fig. [5] we present the ratio R = A 



sm{<t>-<l>s)j^ 
UT 



Dp^ 



/A 



sm{4>-<l>s)j^ 
UT 



This ratio 



ppi 



(which is the same for all three used parametrizations for Sivers function) changes from 0.4 



^^That is why the equaUty f-^j^ 



-L(l)« 



JIT 



in the parametrizations (p8)) - ([30)) is apphed. Together with 



the assumption ([51]) it denotes that the equahty fj^j^ '^ + /^^^i ' = for the sea PDFs is also holds for these 



parametrizations . 

-•^■^Evolution model, which is consistent with the analysis of Ref. |9j (see discussion around Fig. [J), also 
predicts that hiu j^ hid- 
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sm{(f>+4>s)jf-- I -f 

to 0.8. In Fig. [6] we present the asymmetry A^rp ^ for Dp' and pp' collisions and it is 

seen that the respective curves almost merge. Thus, one can conclude that in the case of Dp^ 
collisions both, weighted with sin(0 — (ps) and weighted with sin(0 + (ps), SSA are presumably 
measurable in the same x regions as the respective SSA in the case of pp^ collisions. 

5 Estimations on the SSA feasibility with the new gen- 
erator of polarized DY events 

Generator of polarized Drell-Yan events is necessary for first, estimation of the SSA feasibility 
on the preliminary (theoretical) stage (without details of experimental setup) -just this paper. 
Second, as an input for detector simulation software (for example, GEANT [35] based code) 
on both planning of experimental setup and the data analysis stages. Until recently there was 
no in the free access any generator of Drell-Yan events except for the only PYTHIA generator 
|36j . However, regretfully, in PYTHIA there are only unpolarized Drell-Yan processes and, 
besides, they are implemented in PYTHIA without correct qt and cos 20 dependence, which is 
absolutely necessary to study Boer-Mulders effect. Recently did appear the first generator of 
polarized DY events |37j, [38j where qt and angle dependencies are properly taken into account. 

To simplify implementation of the new possibilities in program and to effectively control all 
calculations we wrote the new generator of polarized DY events (the details will be published 
elsewhere). The scheme of generator is quite simple and very similar to the event generator 
GMC_TRANS [393 which was successfully used by HERMES collaboration for simulation of the 
Sivers effect in semi-inclusive DIS processes |10]. Briefly, the scheme of DY event generation look 
as follows. First, the generator performs the choice of flavor q of annihilating qq pair and choice 
does the given hadron (for example, polarized) contains annihilating quark or, alternatively, 
antiquark of given (chosen) flavor. It is done in accordance with the total unpolarized DY cross- 
sections for each flavor and each alternative choice for given hadron in initial state (annihilating 
quark or antiquark inside). Then, the variables xp and Q^ are selected according to the part 
of unpolarized cross-section (see, for example, Eq. (1) in Ref. [21]) which does not contain 
the angle dependencies. At the next step the polar angle 6 is selected from sin6'(l + cos^ 6) 
distribution in that cross-section. Then, the Gaussian model for fiq{x, kx) is applied and 
after that the transverse momentum of lepton pair q^p is selected from exponential distribution 
exp(— g|./2)/27r. At the next step and 05 angles are selected in accordance with cos 20, 
sin(0 — 05) and sin(0 + (ps) dependencies of the single-polarized DY cross-section (see, for 
example, Eq. (2) in Ref. [21]). The fc-r dependencies of Boer-Mulders h^{x,kT) and Sivers 
flx^x, kx) PDFs are fixed by the Boer model [16] and Gaussian ansatz [71 E[, respectively. At 
this stage of and 05 selection xp, Q^, and qp variables are already fixeco that essentially 
increases the rate of and 05 selection. All variables are generated using the standard von 
Neumann acception- rejection technique (see, for example [36j). 

Let us stress once again that the generator elaborated in [371 EH] is the first generator of 
polarized DY events and in many respects it helped us to construct our generator. Now we 
briefly consider the advantages of the new generator. As it was discussed above, one of the 
main requirement to any generator of polarized DY events is to properly include the nontrivial 
qr dependence of DY cross-sections. This is of especial importance for the qr weighted objects 



^''Certainly, one can select all variables simultaneously. However, such scheme essentially decreases the rate 
of events generation. 
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we deal with in the paper. In the earher generator [371 [38] the approximations are apphed for 
the convolution calculations in the hi containing parts of DY cross-section just as it was done 
in the original paper [16] (see the discussion around Eqs. (47) and (53) in Ref. [16]). However, 
the direct calculations show that implementation of these approximation in the generator leads 

to essential distortion of the qr weighted objects, such as SSA A^j, '^ : the values of SSA 

obtained from simulated data significantly differ from the respective values calculated directly 
from the input parametrizations/models for fig, hi^ and hig. That is why we avoid the such 
kind of approximations. Instead, we calculate the convolutions numerically for the large discrete 
set of qt values and then we perform the standard spline interpolation procedure to reconstruct 
the calculated convolution as the continuous function of qr- As a result (see Figs. 7-10 below), 

the values of SSA A^rp ^ reconstructed from simulated data are in a good agreement 

with the respective values calculated directly from the parametrizations/models entering the 
generator as an input. 

Another advantage of the new generator is that it can be much more easy combinecO with 
the PYTHIA generator (where almost all possible processes are included) which is necessary 
for the analysis of background processes that can produce false DY events (misidentification of 
lepton pair). The point is that, on the contrary to generator [371 [38], iii our generator, just as in 
PYTHIA, the DY processes are generated for each flavor in separation. Besides, the constructed 
generator have some technical advantages: due to the chosen generation scheme the rate of 
event generation is much higher than for generator [STJ [38j , where all kinematical variables are 
thrown simultaneously; the search for cross-sections maximums for the von Neumann algorithm 
is performed automatically as well as the calculation of the total cross-section at the end of 
each generation run. 

Having in our disposal the generator of polarized DY events, we can now estimate the 
feasibility of SSA calculated before. Certainly, these are very preliminary estimations on the 
first (theoretical level). To perform the comprehensive feasibility estimations one needs to take 
into account the all peculiarities of the concrete experimental setup. 

We prepared two samples with applied statistics lOOK and 50K of pure Drell-Yan events 
for each of two Q^ ranges: 2 < Q^ < 8.5 GeV^ and Q^ > llGeV\ Cut 2 < Q^ < 8.5 GeV^ 
is applied to avoid misidentification of lepton pairs due to numerous background processes 
(combinatorial background from Dalitz-decays and gamma conversions, etc - see, for instance, 
section F.4.2 in Ref. ^Tj) below Q^ = 2GeV'^ and to exclude lepton pairs coming from J/ip 
region. Cut Q^ > 11 GeV^ is also applied to avoid the lepton pairs coming from J/ip region. 

As before, we did not present here the estimation on feasibility of Atjj. ^ for RHIC 

ppT 

kinematics because it was in detail done in Ref. pj. So, for this SSA we again (see Section 3) 
present here the results only for the NICA center of mass energy 20 GeV having in mind that 
the results on SSA for COMPASS and J-PARC (in comparison with RHIC rather close in y/s 
to NICA) occur quite similar (see Section 3) to the respective results for NICA. The results 
are presented in Fig. [71 For the simulations with the developed generator we use the latest 
parametrization, fit HI (solid line in Fig. [7]), from the set f[28l) - (!30l) . Looking at Fig. [7] one 
can see that even at relatively low applied statistics 50K pure Drell-Yan events (bottom part of 

Fig. [7]) there are three presumably measurable points for A^rp ^ in the kinematical region 

Xp — XpT > 0, where this SSA is about 4-6%. Moreover, at applied statistics lOOK pure Drell- 



15 



This work now in progress. 
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sm((f)—cf>s)j^f— 

Yan events one can hope to reconstruct the functional form of SSA Ajjrp ^ (see top part 

of Fig. [7]) in the kinematical region Xp > XpU Since kinematical region Xp > Xpt corresponds 
to xj? > for the fixed target mode with unpolarized beam and polarized target available to 
COMPASS and J-PARC, one can conclude that all four RHIC, NICA, COMPASS and J-PARC 
facilities can provide us the access to the Sivers PDF (see Section 3). In the region Xp < Xpt 

sin(0— (/)g) j^j^ 

available to RHIC, NICA and (presumably) J-PARC, SSA A^J^ '^ is smaller (less than 

4%), but still visible within the errors (even at applied statistics 50K events one can see at least 
one measurable point). 

Let us now consider feasibility of the Sivers PDFs in the case of the deutron in initial state. 

s 

Looking at Figs. [5]and[7]one can conclude that the only SSA Ajjj. 



.sin(0-</,s)^ 



which could take 

Dpt 

considerable values (see Section 4) is hardly measurable at applied statistics 50K pure Drell- 
Yan events. At the same time, in the case of lOOK events this asymmetry becomes presumably 
measurable. 

sm{<j)+<t,s)jf- 

Let us now estimate the feasibility of SSA A^rp ^ giving us an access to transversity 

and Boer-Mulders PDFs. We again (see Section 3) present the estimations for pp^ collisions 
and for two quite different RHIC and NICA center of mass energies. The results are presented 
in Figs. [H]and[ni For COMPASS and J-PARC kinematics our calculations produce almost the 
same figures as for NICA center of mass energy, so that we again omit the respective plots. For 
the simulations with the developed generator we again use (see Section 3) Boer model for hj^ 
and the evolution model for hi with /iig(g) = Ag(Ag) ansatz at initial scale Ql = 0.23 GeV^. 

Looking at Figs. M and M one can see that in the region Xp < Xpt even at statistics 50K 
pure Drell-Yan events (bottom part of Fig. [H]) one can hope to see within the errors at least 

sm{<P+<l,s)jf- 

three points for A^rp ^ in the kinematical region Xp — Xp-\ < 0. At the same time, at the 

statistics lOOK events one can hope also to reconstruct the functional form of this SSA (see 
top part of Fig. [8]) in the kinematical region Xp < Xp^. Regretfully, the kinematical region 
Xp < Xpf is hardly accessible for COMPASS because of the forward geometry spectrometer 
and the unpolarized proton beam available to this running experiment. In the case of another 
experiment with the fixed target, planning at J-PARC facility, one could reach this kinematical 
region if the polarized proton beam would be available (this option is now planned at J-PARC 
[2]). Fortunately, there are no any problems to reach the kinematical region Xp < Xpt for the 
symmetric collider mode available to RHIC and NICA. 

Let us note that due to the close values of Ajjj. '^ in the cases of pp^ and Dp^ collisions 

sin((l)+4>s)jf- 

(see Fig. E]) all conclusions concerning feasibility of Ajjj, ^ in the case of pp^ collisions 

are valid in Dp^ case too. 

Looking at Figs. [319] one can see that if the available experimental statistics will be rather 
low (less than 50K events) it is hardly possible to reconstruct the functional form of Sivers and 
Boer-Mulders PDFs. However, even in this unfavorable case one can hope at least to check 
very important QCD prediction ([T]). 

siyi{<f>—4>s) tT sin((/)+(i!)g) ^ 

In conclusion of this section let us estimate the feasibility of both A^j, ^ and Ajjr^ ^ 

SSA for the n~p^ —>■ /i+/i~X DY processes available to COMPASS. The model estimation of 

SSA A^rp '^ magnitude was already presented in the previous paper f^2\- However, that 

time we could not present the feasibility estimations for this SSA because there still was not 
any generator of polarized DY events in the free access. Besides, only recently the optimal 
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kinematical conditions for DY program at COMPASS were chosen. These are the pion beam 
energy 160 GeV and the range 4 < Q < 9 GeV for the invariant dilepton mass [3] . The results 

on SSA A^rp ^ _ are presented in Fig. [101 Because of the restriction f[32|) for the 

COMPASS forward geometry spectrometer we present the results only for x.^ > XpU Looking 
at Fig. [in] one can see that even at rather low applied statistics 50K events there six presum- 
ably measurable points for both asymmetries, which encourage us that it could be possible to 
reconstruct the functional form for both SSA. Thus, one can conclude that for COMPASS it 
is much more profitable to study DY processes with the pion-proton collisions. At the same 
time, as it was discussed above, DY processes p{D)p\D^) — ^ l^l~X could be most successfully 
studied in collider mode at RHIC and NICA. 

6 Conclusion 

In summary, the Drell-Yan processes with the colliding protons and deutrons available to RHIC, 
NICA, COMPASS and J-PARC were considered. We estimated the single-spin asymmetries 

A^J^ ^ and A^J• ^ , which give us an access to Sivers and to Boer-Mulders and 

transversity PDFs, respectively. The estimations were performed for the different y/s values, 
corresponding to the kinematical conditions of RHIC, NICA, J-PARC and COMPASS facilities. 

The preliminary estimations for pp' collisions demonstrate that SSA Ajjrp ^ can reach 

quite considerable values (5-10%) in both Xp > Xp] and Xp < Xpt regions. On other hand, 

the estimations performed for SSA A^J• ^ show that this asymmetry is negligible in 

the region Xp > Xp^ and takes considerable values (also about 5-10%) in the region Xp < 
XpT. The estimations performed for SSA in the case of Dp and DD collisions demonstrate 
that the asymmetries for DY processes with pD^ collisions are compatible with zero except 
for perhaps one, containing sum hiu + hid- The latest analysis [9] predicts that this sum 
could essentially differ from zero. Certainly, because of very large uncertainties, this is very 
preliminary conclusion which should be carefully checked in the future DY experiments. On 
the contrary to DY processes with pD^ and DD^ collisions, SSA for Dp^ collisions are close in 
their values to the respective SSA for pp^ collisions and, thus, presumably could be feasible in 
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the same kinematical regions. While SSA Ajjrp 



the asymmetry Arjrp ^ is even slightly larger than Arjrp 

DpT ^ pp\ 

The new generator of polarized DY events was developed which allowed us to estimate the 
feasibility of both weighted with sin(0 — (ps) and sin(0 -|- (ps) single-spin asymmetries. These 
estimations performed for proton-proton collisions demonstrate that both SSA are presumably 

measurable even at the applied statistics 50K pure Drell-Yan events. While A^j, ^ is 

presumably measurable in both kinematical regions Xp > x^t and Xp < x^t, the asymmetry 

sm{(f>+4>s)jf- 

Auj, ^ could be measured only in the region Xp < XpT, where it takes quite considerable 

values. It is of importance that while for the symmetric collider mode (RHIC, NICA) there 
is no any difficulties to reach both Xp > Xpt and Xp < x^t kinematical regions, for the fixed 
target mode (COMPASS, J-PARC), there is the problem to reach the region Xp < XpU This is 
hardly possible for the running COMPASS experiment, where only unpolarized proton beam 
and forward geometry spectrometer is available. At the same time it still could be achieved by 
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J-PARC, where at present the option with polarized proton beam is planned. 

We also studied the behavior of SSA in two different limiting cases, Xp ^ Xp] and Xp -C 
XpT . It is of importance that in these limiting cases one can essentially reduce the number of 
unknown PDFs entering the asymmetries. In particular, studying the unpolarized and single- 
polarized Drell-Yan processes in the limiting case Xp <^ x^t one can directly extract the ratio 
of transversity and first moment of Boer-Mulders PDF. 

Thus, one can conclude that it is much better to study the single-polarized DY processes 
with pp, pD and DD collisions in the symmetric collider mode (RHIC and NICA). For fixed 
target experiments (like J-PARC and COMPASS) the polarized proton beam is necessary to 
improve the situation. This option is planned at J-PARC but is hardly possible for already 
running COMPASS experiment, where it is much more profitable to study DY processes with 
the pion-proton collisions (namely this option now is planned at COMPASS |3]). 

Let us stress once again that any new measurements of SSA in Drell-Yan processes are 
of extreme importance even at rather poor available statistics of collected Drell-Yan events. 
Even in this unfavorable case one can at least to check very important QCD prediction ([T]), 
which would provide a crucial test of our understanding of T-odd effects within QCD and the 
factorization approach to the processes sensitive to transverse parton momenta. It is encourag- 
ing that today the Drell-Yan measurements are planned simultaneously in the different world 
centers for high energy physics. 

Acknowledgments 

The authors are grateful to R. Bertini, O. Denisov, A. Efremov, Y. Goto, T. Iwata, S. Kazutaka, 
V. Kekelidze, S. Kumano, A. Maggiora, I. Meshkov, A. Olshevsky, G. Piragino, G. Pontecorvo, 
S. Sawada, I. Savin, A. Sorin and O. Teryaev, for fruitful discussions. 

The work of O. Shevchenko, A. Nagaytsev and O. Ivanov was supported by the Russian 
Foundation for Basic Research (Project No. 07-02-01046). 

References 

[1] RHIC Spin Collab., D. Hill et al., letter of intent RHIC-SPIN-LOI-1991, updated 1993; 
G. Bunce et al. Particle World 3, 1 (1992); 

PHENIX/Spin Collaboration, K. Imai et al., BNL-PR0P0SAL-R5-ADD (1994). 
O. Martin, A. Schafer, M. Stratmann, W. Vogelsang, Phys. Rev. D57 (1998) 3084; Phys. 
Rev. D60 (1999) 117502. 

[2] J. Chiba et al, J-PARC proposal "Measurement of high- mass dimuon produc- 
tion at the 50-GeV proton synchrotron", can be obtained electronically via 
http://j-parc.jp/NuclPart/pac_0606/pdf/p04-Peng.pdfi 



[3] M. Colantoni, talk given at Transversity 2008 workshop. May 28th - May 31st, 2008 
Ferrara, Italy; accessible electronically via http://www.fe.infn.it/transversity2008/| 

[4] A.Sissakian et al.. Conceptual Design Report "Design and Construction of Nuclotron- 
based Ion Collider facility (NICA)", Dubna 2007. Accessible electronically via 
[http://nucloserv.jinr.ru/nica_webpage/Nica_files/reports/CDR_07/CDR_NICA^ 20.html 

[5] A.N. Sissakian, A.S. Sorin and V.D. Toneev, QCD Matter: A search for a mixed quark- 
hadron phase, |nucl-th/0608032[ 

16 



[6] V. Barone, A. Drago, and P.G. Ratcliffe, Phys. Rep. 359, 1 (2002). 

[7] A. V. Efremov et al, Phys. Lett. B612 (2005) 233 

[8] J.C. Collins et al, Phys. Rev. D73 (2006) 014021 

[9] M. Anselmino et al, Phys. Rev. D75 (2007) 054032 

[10] HERMES collaboration (A. Airapetian et al), Phys. Rev. Lett. 84, 4047 (2000); Phys. 
Rev. D 64, 097101 (2001); Phys. Lett. B 562, 182 (2003); Phys. Rev. Lett. 94, 012002 
(2005). 

[11] E.S. Ageev et al (COMPASS collaboration), Nucl. Phys. B 765 (2007) 31 

[12] J. C. Collins, Phys. Lett. B 536 (2002) 43 

[13] J. Collins, Rapidity divergences and valid definitions of parton densities. Proceedings of 
LIGHT CONE 2008; eprint: .arXiv:0808.2665v 2 [hep-ph] 

[14] V. Barone, A. Prokudin, B.-Q. Ma, Phys. Rev. D78 (2008) 045022 

[15] O. Martin, A. Schafer, M. Stratmann and W. Vogelsang, Phys. Rev. D57 (1998) 3084 

[16] D. Boer, Phys. Rev. D 60, 014012 (1999) 

[17] A. Kotzinian, P.J. Mulders, Phys. Letts. B 406 (1997) 373. 

[18] D. Boer, R. Jakob, P. J. Mulders, Nucl. Phys. B 504, 345 (1997) 

[19] D. Boer, R. Jakob, P. J. Mulders, Phys. Lett. B 424, 143 (1998) 

[20] D. Boer, P. J. Mulders, Phys. Rev. D 57, 5780 (1998) 

[21] A.N. Sissakian, O.Yu. Shevchenko, A.P. Nagaytsev, O.N. Ivanov, Phys. Rev. D72 (2005) 
054027 

[22] A. Sissakian, O. Shevchenko, A. Nagaytsev, O. Denisov, O. Ivanov Eur. Phys. J. C46 
(2006) 147 

[23] J.C. Collins et al, Phys. Rev. D73 (2006) 094023 

[24] M. Cluck, E. Reya, A. Vogt, Z. Phys. C67 (1995) 433 

[25] M. Anselmino et al, Phys. Rev. D 72 (2005) 094007 

[26] W. Vogelsang, F. Yuan, Phys. Rev. D 72 (2005) 054028 



[27] A. Bacchetta, U. D'Alesio, M. Diehl, C. A. Miller, |arXiv:hep-ph/0410050 

[28] NAIO Collab., Z. Phys. C31, 513 (1986); 

[29] J.S. Conway et al, Phys. Rev. D39, 92 (1989). 



[30] V. Y. Alexakhin, et al., Phys. Rev. Lett. 94, 202002 (2005), hep-ex/0503002 



17 



[31 
[32 
[33 
[34 

[35 
[36 

[37; 

[38 
[39 

[40 

[41 
[42 
[43; 



Belle Collaboration, R. Seidl et al, Pliys. Rev. Lett. 96, 232002 (2006). 

V. Barone , T. Calarco, A. Drago, Phys. Rev. D56 (1997) 527 

P.V. Pobylitsa, h ep-ph/0301236] 

A. V. Efremov, K. Goeke and P. V. Pobylitsa, Phys. Lett. B 488 (2000) 182 
[a rXiv:hep-ph/0004196j . 



See http://wwwasd.web.cern.ch/wwwasd/geant/ 



T. Sjostrand et al., hep-ph/0308153 



A. Bianconi, M. Radici, Phys. Rev. D73 (2006) 034018; Phys. Rev. D72 (2005) 074013 



A. Bianconi, arXiv:0 806.0946 [hep-ex] 

N.C.R. Makins, GMCJrans manual, HERMES internal report 2003, HERMES-03-060; 
G. Schnell, talk at workshop Transversity'07, ECT* Trento, Italy, June 2007 

A. Airapetian et al. (HERMES Collaboration), Phys. Rev. Lett. 84, 4047 (2000); Phys. 
Rev. D 64, 097101 (2001); Phys. Lett. B 562, 182 (2003) 



V. Barone et. al (PAX collaboration), hep-ex/0505054 



M. Cluck, E. Reya, M. Stratmann, W. Vogelsang, Phys. Rev. D53 (1996) 4775 
M. Cluck, E. Reya, A. Vogt, Z. Phys. C53 (1992) 651 



sin(i^-<^, 




sin(<^-y^ 



0.14 


r 




11 




^^^ 


0.12 


r 





III 




/ 


0.1 


r 








/ 

/ 


0.08 


r 








//;^j~^^^^^ ^ 


0.06 
0.04 
0.02 


r 




^ 


y"^ 








rrr 


-jTp-'^ 


1 , , 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



-0.8-0.6-0.4-0.2 -0 0.2 0.4 0.6 0.8 

Xp-XpT 



Figure 1: Estimation of SSA A^t ~ ' ^\ ^ for NICA, s=AOOGeV\ with Q^ = AGeV^ (left) 

and Q^ = 15GeV'^ (right). Rome numbers /, // denote respectively fits I and II from Ref. [^ 
and III denotes the fit from Ref. [S]. 
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Figure 2: Results of Ref. p] on hi^ in comparison with the resuhs obtained with two versions 
of evolution model and with the Soffer bound. The bold solid line (top) corresponds to upper 
bound given by the Soffer inequality. Dashed line corresponds to the evolution model with the 
Soffer inequality saturation at the initial model scale Qo = 0.23 Gel^^. The solid line corre- 
sponds to the upper boundary of the error band on hiu- The dashed-dotted hne corresponds to 
the evolution model, where hiu^u = Au{Au) at initial scale Ql = 0.23 GeV^. The dotted line 
corresponds to the fit of Ref. [9] on hiu- GRV94 [21] parametrization for q{x) and GRSV95 
parametrization for Aq{x) are used. 
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Figure 3: Estimation of SSA A 
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(left) and Q^ = 20GeV^ (right). The solid and dotted curves correspond to the two different 
input ansatzes for hi^ which are used in evolution model. These are hig^g = Aq, q and hig = 
{Aq + q)/2 hiq = {Aq + q)/2, respectively. Here GRV94 [21] parametrization for q{x) and 
GRSV95 [12] parametrization for Aq[x) are used. 
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Figure 4: Estimation of SSA Ajjj, ^ 
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for NICA, s = AOOGeV^, with Q^ = AGeV^ (left) 



and Q^ = IbGeV"^ (right). The sohd and dotted curves correspond to the two different input 
ansatzes for hiu which are used in evolution model. These are hig^g = Aq, q and hig = {Aq+q)/2 
hiq = (Ag + q)/2, respectively. Here GRV94 [21] parametrization for q{x) and GRSV95 
parametrization for Aq{x) are used. 
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Figure 5: Estimation of ratio R = A 

Q2 = 4GeV^ (left) and Q^ = i^GeV^ (right). 
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Figure 6: Estimation of SSA A^jj. '^ for NICA with two different options, pp^ (dashed 

line) and Dp^ (sohd line) colhsions, with Q"^ = AGeV"^ (left) and Q"^ = 15GeV'^ (right). Here 
evolution model for hig is used with the input ansatz /iig g = Aq, q. GRV94 [2^ for q and 
GRSV95 [l2] parametrizations for Ag are used. 
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Figure 7: Estimation of asymmetries A^rp 

Ref. [H] is used. The points with error bars are obtained by using simulations with event 
generator at the apphed statistics lOOK pure Drell-Yan events. (Q^) ^ 3.5 GeV^ and {Q"^) 
15 GeV^ for the left and right plots, respectively. 
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Figure 8: Estimation of asymmetries Ajjj^ ^ 

evolution model with the input ansatz hiq^q = Aq,q at Ql = 0.23 GeV^ is used. GRV94 
[M] parametrization for q{x) and GRSV95 |12] parametrization for Aq{x) are used. The points 
with error bars are obtained with the developed generator of polarized DY events at the applied 
statistics lOOK (top) and 50K (bottom) pure Drell-Yan events. {Q"^) ~ 3.5 GeV'^ and {Q"^) = 
15 GeV^ for the left and right plots, respectively. 
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Figure 9: Estimation of asymmetries Ajjj, ^ 

evolution model with the input ansatz hiq^q = Aq, q at Ql = 0.23 GeV^ is used. GRV94 
|24j parametrization for q{x) and GRSV95 |12] parametrization for Aq{x) are used. The points 
with error bars are obtained with the developed generator of polarized DY events at the applied 
statistics lOOK (top) and 50K (bottom) pure Drell-Yan events. (Q^) ~ 3.9 GeV^ and {Q^) = 
22 GeV^ for the left and right plots, respectively. 
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Figure 10: Estimations on SSA A^rp 
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